Symmetry in Order-Disorder Changes of Molecular Clusters 



Ana Proykova^, Dessislava Nikolova^, and R.Stephen Berry* 
^University of Sofia, Faculty of Physics, 5 James Bourchier Blvd., Sofia-1126, Bulgaria 
* The University of Chicago, Department of Chemistry Chicago, IL 60637, USA 

The dynamic orientational order-disorder transition of clusters consisting of octahedral AFa 
molecules is formulated in terms of symmetry-adapted rotator functions. The transition from a 
higher-temperature body-centered-cubic phase whose molecules are orientationally disordered at 
their sites to lower-temperature, monoclinic, orientationally-ordered phase is a two-step process: 
first, at temperatures well below the limit of stability for the liquid, a transition occurs to a partially 
ordered monoclinic phase driven by the rotational- vibrational coupling. This transition has two local 
minima in the free energy, and hence behaves like a finite-system counterpart of a first-order transi- 
tion. Further lowering of the temperature initiates another transition, to an orientationally-ordered 
base-centered monoclinic structure. This last transition is dominated by rotational-rotational inter- 
action and is found from simulations to be continuous. The temperature of this transition predicted 
by the analytic theory presented here for a 59-molecule cluster of TeFg, 27K, is in good agreement 
with the 30K result of canonical Monte Carlo calculations. 
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Below the freezing point, the plastic phases of molecular substances are known to form highly symmetrical lattices 
with orientationally disordered molecules ^ . As the temperature of the disordered substance is still lowered, a phase 
transition to a more orientationally-ordered phase takes place. The new phase might have a crystal structure of lower 
symmetry, characterized by long-range orientational order. This is called displacive-ordering transition. Similar struc- 
tural transformations have been also detected in small, free clusters consisting of rigid octahedral molecules SFq Q| 
and TcFq and have been studied by numerical simulations and experimentally . The change between 

phases may be continuous or discontinuous depending on the interaction potential. However, simulations of clusters 
always produce smooth behavior in the sense that clusters in different phases coexist over a range of temperatures 
and pressures, regardless of the order of the transition in the thermodynamic, large-N limit. In many cases, the 
crystal forms observed depend strongly on the conditions of production, thus leading to different interpretations of 
what should be the final structures of the clusters at low temperatures. Structures such as triclinic, monoclinic and 
hexagonal were reported for the same material In our previous simulations, we have shown that the finite size 
of clusters causes rotational [Q and strong surface effects |lj], appearance of vibrational-rotational coupling, none of 
which are expected in bulk systems of the same symmetry. 

Vibrational-rotational coupling was considered qualitatively to explain the two-step process of ordering in some 
clusters of octahedral molecules, Oh symmetry, as the temperature decreases [|[^. Below the melting temperature, 
those clusters assume a body-centered cubic structure bcc, Oh symmetry, and orientational disorder of the molecules 
Simulations show that this transition involves passage between two forms at different local minima in 
the free energy, and hence is the apparent counterpart of a first-order transition. Whether these two minima remain 
apart or converge as A'^ — > oo is an unsettled question [^,^ . A temperature decrease drives a phase transition from the 
orientationally disordered bcc structure to an orientationally ordered monoclinic phase, C2/1. Since no normal mode of 
Oh carries the group directly to C2h this transition could occur in two steps: first, the bcc transforms into a partially 
ordered monoclinic phase driven by rotational-translational coupling. Second, after a further temperature decrease, 
another transition to an orientationally ordered phase without change of the symmetry of the lattice sites takes place. 

In the present paper we analyze analytically the dynamic orientational order in clusters of rigid octahedral molecules 
of type AFq to determine the contribution of various interactions to the total potential and the nature of the 
temperature-driven solid-solid transformation of plastic clusters. The symmetry of the molecules, Oh, and the 
molecular sites in the cluster. Oh or C2h, are taken into account explicitly in our analysis, based on the theory 
of orientationally-disordered crystals . The site position of each molecule is determined by the Cartesian coor- 

dinates of its mass center at A. Throughout the paper we refer to (i) a space (laboratory) axis system XYZ, fixed at 
the cluster's center of mass; (ii) a nonrotating system X'Y'Z' parallel to XYZ but with its origin translating with 
the molecular center of mass; (Hi) a rotating, body-fixed system xyz. For rigid molecules, these axes coincide with 
the principle axes of inertia. 

The translation of a molecule is separable as the motion of the molecular center of mass in the XYZ coordinates. 
The xyz orientation in the XYZ system is given with fl = (0, 4>, ip), where 9 and (f) are the ordinary polar coordinates 
of the z axis in the XYZ system and ■0 is an angle in the xy-plane measuring the rotation clockwise about the z axis. 

The molecular orientation in orientationally-disordered high- T crystals has been presented in terms of symmetry- 
adapted rotator functions used by James and Keenan to describe the orientational phase of methane with tetra- 
hedral symmetry. 

Let us consider an octahedral molecule n in its initial orientation Qn where the molecular axes coincide with the 
laboratory system axes. The orientational density distribution is expanded in terms of spherical harmonics YJ™(r2). 
The molecular symmetry requires that I = 0,4,6,.... and only certain linear combinations of YJ™ occur. For each 
allowed /, we determine the molecular symmetry-adapted functions 



where the superscript A refers to the identity representation of the cubic group Oh and a™ are tabulated in Ref. 
A molecular form factor is defined for the allowed ^=0,4,6,...as 



if the molecular axes coincide with the space axes; Na is the number of the atoms in a molecule; Qi,{n) denotes the 
orientation of atom v in the space-fixed system. 




(0.1) 




(0.2) 
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Unlike bulk crystals, the cubic structure of a free cluster is broken at the cluster surface. In what follows, we 
consider the symmetry of the volume molecules, defined as molecules having all the neighbors required by a specific 
point group. We neglect the symmetry-breaking associated with the surface molecules, which belong to a different 
group, if any. The separate description of the surface and the volume limits the analysis. In the closing remarks, we 
discuss a possible correction to this assumption. 

The orientational changes at a cubic site are also expressed in terms of site-symmetry-adapted functions p^ ]: 

I 

s[{n) = Yl Yri^)^r, (0.3) 

m— — l 

where the superscript t = (G, F, p, p) indicates the irreducible representations F of the group G, p distinguishes 
between the representations that occur more than once, and p denotes the rows of a given representation. The 
symmetry-adapted functions represent a complete basis in the ^{9, (/), 7/')-space. The I = 4 manifold of this system 
reduces to the representation Aig, Eg, Tig, T2g under G = Oh and to the representations bAg and ABg under G = G2/1. 
For the normalized function S'[ we find: 

= [0.763, m = 0; 0.457, m = ±4] for Oh, Aig, 1, 1 
xi [0.645, m==0; -0.541, TO =±4] for Eg,\,l 
af^ = [-0.707, TO = ±2] for Oh, Eg, 2, 1 
af^ = [-iO.663, m = ±1; -iO.25, m = ±3] for Oh, Tig, 1, 1 
= [±i0.663,m= ±l;-i0.25,m= ±3] for Oh, Tig, 2,1 
af^ = [±i0.707, m = ±4] for Oh, Tig, 3, 1 
af^ = [±i0.663, m = ±3; ±i0.25, m = ±1] for Oh, ^29, 1, 1 
af^ = [±0.663, TO = ±3; ±0.25, TO = ±1] for Oh, ^2^,2,1 
af = [±i0.707, m = ±2] for Oh, T2g, 3, 1. 
The largest contribution to the crystal field is that of the Aig component of the I = 4 manifold. 

At low temperatures, clusters of AFq molecules adopt a monoclinic structure C2h- AH coefficients are equal to 
1 for the five-fold representation Ag {m = 0, ±2, ±4) and for the four-fold representation Bg (m = ±1, ±3). 

An arbitrary molecular orientation f2„/ with respect to the initial one f2„ is obtained by a rotation specified with 
the Eulcr angles (w=a,/?, 7). Rotation fj'' of a molecule does not affect the spherical harmonics Y;'=(f7'') defined in the 
body system. In the space system these are changed to R{lo)Y^{U!') = X]L=-/ Yi^{^'^)DY^'^{lo), where SI* determines 
the space orientation of the molecule. For a molecule in arbitrary orientation f2„' , the symmetry-adapted function 
changes to 

k— — l 7n— — l 
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with D]^ (lu) the Wigner matrices. We determine the spherical harmonics Fj™(ri'') from the equation inverse to (0.1) 
and put it in ( |0.4| ). The result is 

R{co)s^in'') = ^ Sim (aTl DrH^hf' (0.5) 

The equation ( |o7^ ) relates the symmetry-adapted functions S^{Q'') for the body (molecular) system and the 
symmetry-adapted functions Si{i}^) for the space (cluster) system. Rewriting the cq. ( |0.5D as R{u!)Sl'{i}'^) = 
Y^^Sl{n'')AJ^{u;), where A[^(a;) = J2k=-iT,L=-i (a™^) are the rotator functions A[^(w) defined by 
the symmetry properties of the molecule af^ and of the site a™"^ . Rotator functions with I = 3 were introduced for 
solid CD4 in [^. The rotator function's average value A.J-^{ll!) is zero in the disordered phase and non-zero in the 
orientationally ordered phase. This property makes it suitable to be chosen as an order parameter. 

So far we have considered a single molecule at a specific site. The orientational configuration of N molecules in the 
cluster is given by AJ''^(w(n)) = AJ^{n) where n = 1, 2, N, labels each molecule's center at its lattice position r„. 
The interaction between two molecules n and n' can be written as a sum of atom-atom potentials ^ : 
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where (n, i/) labels the ly*^ atom in the molecule at site r„. The total potential V oi N molecules is: 

N Na 



n<n' i/.y' 

The potential V{n, v] n'v') depends on the distance r^^i between the atoms v and v' . The position of the z^*'' atom 
in the n*"^ molecule with respect to the space system is given by 

'R{n,(F) = r„ + (fn^{n) + u(n) 

with u(n) being the displacement of the n*-^ molecule from its site position r„. Vl^{n) indicates the orientation of the 
vector r„y in the space system and (F is its length. We expand V from Eq.(p.6[) in terms of the displacements u(n): 



V{n,n') = EE 7^^n''-p(^-') ["'^ ^ " - " (0.8) 

p=0 v,v' 

with the notation 

The coefficients contain the orientational dependence of the molecules at the sites n and n' . We expand them 
in terms of symmetry- adapted functions SJ ( |0.3| ). In the following we write S'^(j^) for SJ {Vl^), where = //(r, I): 

mm' 



^n...v('^-') - E^if-.MM'(^'"')^MM5,'(^')- (0.10) 



The coefficients cf^ ■ ,.,,An,n') are determined from 



ii...ipfifi' 



,{n,n') ^ dn, dn,.Vl^^^^{n,n')S^,{u)S^,{v'). (0.11) 



We put V^P^ from Eq.( |0.10| ) in Eq.(l^) and use the molecular form factor Eq.( |0.2| ) if the molecular axes coincide 
with the space axes or X]i/'^m('^) ~ di^^tJ.i'^) if the molecule is rotated at an angle uj. Thus the pair potential becomes: 

V{n,n') = EE ^<.-pMM'('^'"')3^i'A,,(n)A,,,(n') [uM ~ n^M)] - - u,^(n')] (0.12) 

p mm' 

Molecular and site symmetry considerations restrict the number of terms in the sums, thus redu cing t he computa- 
tional effort to obtain the contribution of the different interactions in the total pair potential Eq. ( 3.12 ). 

The value of p = corresponds to a rigid lattice (no displacements of the molecular center of mass) . For this case 
we get only rotational-rotational interaction between two molecules with /i ^ (0, 0) and ^' ^ (0, 0): 

l/"(n,n') = Ec°,-(","').9'.9rA4n)A^,(n') (0.13) 
mm' 

The total rotational interaction is the sum over all molecules Ylinn' ^"("j "n')- 
The matrix of rotational-rotational interaction is defined by 

J^^\n,n')^cZ{n,n)gUv (O-W) 

where 

cf^'iri^n) ^ Jdn^J dn,,vl'^\{n,n')SMS^.'{i^')- (0.15) 

The structure of the interaction matrices cjj'j], (n, n') depends on the symmetry of 5^ and on the relative position 
(n, n') of two interacting molecules on a lattice with a symmetry specified by r. 
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The cluster transforms from a disordered cubic to an ordered monoclinic structure at Tc that is the temperature 
at which the free energies of the two forms are equal. In order to calculate Tc we also need the total field acting 
on the molecule at site r(n). The zero*'* approximation is to consider spherical-symmetrical molecules /z' = (0,0) 
acting on a molecule n ^ (0,0) on a rigid lattice p = 0: V^(°)(n,n') |;'=o= ^^i^ii^ Setting 
p! = (0,0) yields 3^,,= Sq'" = (47r)-i/2 and go = Na{4n)-^/^. The coefficients c^°\n,n') become: c';°\n,n') = 

Let us denote the interaction matrices cj^\n,n') weighted with the molecular factors and gQ with v^, where a 



is an index for {I, Aig,p, p): v^^ = J2n' '^^ ("-i "Off; 5o •• The crystal field acting on the molecule n is: 



(0.16) 



The rotator functions (uj) in Eg. ( 0.16] ) are cubic functions. 

The pair vibrational-rotational interaction is obtained from (PI ) forp= 1, ^ 7^ (0,0), and p' = (0,0): 



y™(n,n') = 1/(1) (n,n') = ^ c^;^) (n, n')5;'5o'A^(n) [u,{n) ~ u,{n')] . 



Ifl 



(0.17) 



The sum over all molecules results in the total bilinear interaction V^'^ = ^„ V^^{n,n'). 

The translational-orientational interaction is caused by the change of the orientational potential due to the dis- 
placement of the nearest neighbors. 

The pair vibrational-vibrational interaction is obtained from Eq.( |0.12 ) for p ~ 2, p' ~ (0,0), p — (0,0): 



V^^{n,n') = V^^\n,n') = ^ \cfX{n,n')g^g^, [uM " u.M] K (n) - u^] 



(0.18) 



This gives for the total vibrational- vibrational interaction V'^'^ ~ ^^^^ V'^'^{n, n'). Now the total potential is: 



RR 



(0.19) 



The equation ( 3.19 ) may be expanded with higher order terms which may become important in some structural phase 
transitions ||l7[| . 

Having determined the interactions and the total field, we can calculate the free energy F of each phase as a function 
of the rotator functions considered as order parameters [p^ : 

F — 0.5 ixo^ + FT[J] d^{q)d^'{~q), where FT(J) and ^^(g) are the Fourier images of the rotator matrix J 

and A^(w), respectively; i is the 3x3 unit matrix; xo = xT~^ is the single molecule orientational susceptibility [ p^ : 
X — Z~^gi J dhJexp{—V^'{uj)/T) (A^(a;))^ with Z = j dujexp{—V^/T) the partition function. The expectation value 
of X does not depend on the components of the rotator function A. Two phases of clusters coexist in equal amounts 
or with equal frequency when their free energies are equal. In the limit of — > cxd, a phase transition occurs at Tc, 
which is the point where an eigenvalue of [IT -|- x{T)¥'l[J\ vanishes, if the transition is continuous. The temperature 
dependence of x is very weak which means that the Curie- Weiss law xo — x{T)/{T — Tc) is valid for negative diagonal 
elements of J. The transition point Tc occurs at the largest value of the matrix for the representations allowed by the 
symmetry of the system, i.e. Tc = max[—xJ]. 

As an example, we have applied this group-theoretical approach to the case of a cluster containing 59 TeTg molecules 
and have compared the results to those published in . In order to account for the broken symmetry at the cluster 
surfaces, we use sum over the nearest neighbors. Since the molecules of type AFq have no low-order electrostatic 
moments, the Coulomb contribution to the Lennard- Jone s pot ential can be neglected, see fig. 2 in Q. 

We compare the rotational-rotational contributions, Eq.( 0.13| ), with their vibrational-rotational counterparts, (D.17), 
for the TcFq molecule with the nearest neighbors located at sites with cubic or monoclinic symmetry. For the purpose, 
we determine the vibrational spectrum, the crystal field ( 0.16| ), and the rotational matrix ( 0.14 ). The CERN Library 
pO| is used to compute he elements of ( 0.14 ). The contribution from I — I' = 4 is the largest. The J matrices in a 
cubic symmetry environment is: 



/ 0.011 -0.0007 \ 
JT2g = -.0007 -0.02 

\ 1.63/ 
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/ -0.004 -0.001 \ 
Jrig = -.001 -0.011 

V -3.81 / 



0.17 -0.02 
-0.02 2.41 



These matrices are diagonal in C2h because it has only one-dimensional representations. The diagonal elements of 
J in the case of C2h are: 

J(Xn) = [0.043, 0.043, -0.027, -0.027, 0.014, 0.014, -0.011, -0.011] 
^(^12) = [-0.41, -0.41, 0.018, 0.018, -0.003, -0.003, -0.002, -0.002] 
^(^13) = [-0.22, -0.22, -0.08, -0.08, -0.0005, -0.0005, -0.0004, -0.0004] 

The largest value of J is obtained for the representation Aig, so we choose as an order parameter the rotator 
functions for this representation and expand the free energy F in terms of the corresponding rotator functions. Aig 
is the only common representation for Oh and C2/1. To find correlation of Oh and C2h we carry out the correlation in 
two steps: first, pass from Oh to D^h and then imply the table for D^h to go on to C2/1 [ |l5|]. 

The coefficients necessary to compute the crystal field (D.16) are calculated from Eq.( p.l5 ). We get for the cubic 

Oh 



symmetry Vq — —7.98 in the approximation of the nearest eight neighbors. For the monoclinic structure, C2h, this 



coefficient is Vq — —0.61. From Eq.(3.2) we obtain for the octahedral molecule ofTeFQ-. (7o=l-98 and 174=1.29. 



In the approximation of the nearest neighbors interaction we have determined the energy per molecule: in Oh , the 
rotational-rotational energy is AmeV and the vibrational-rotational is ImeV; in C2/1, these values are l.2meV and 
0.02mey, respectively. The conclusion is that the vibrational-rotational interaction can be neglected in the ordering 
of molecules on monoclinic sites, so that the lower-temperature transition is entirely driven by rotational ordering 
and that transition is continuous even in small systems. However on cubic sites, motion of the molecular centers of 
mass must be taken into account, and the transition acts at least in small systems like it is first-order. The highly 
degenerate state of an octahedral molecule in the octahedral environment is resolved by a distortion of the cluster if 
the model requires a rigid molecule. In this case we must deal with the Jahn- Teller effect that distorts the cluster. 
This implies that a transition to a lower-symmetry structure is initiated. The larger value of the rotational-rotational 
interaction, however, leads to a partial ordering of the molecules which we clearly see [^ in the thermodynamically 
less favored state of clusters of small size. The transition from a cubic to monoclinic structure is resolved with the 
appearance of the representation of Eg in Dh which is equivalent to conden satio n of an active mode. 

Taking into account the rotational- rotational interaction in the total field ( |0.19 ) for the orientational ordering in C2h 



we obtain Tc — max[~xJ] ~ 27 K. This is in a good agreement with the result (~ 3 0^) c omputed in the molecular 
dynamics simulations Q]. Thus we show that the choice of cubic rotator fimctions (3.16) as an order parameter is 
suitable for describing the mechanism of phase changes. 

A final comment is that the theory of discrete point groups used in the present work sets limits on handling at 
once the surface and volume symmetries. However, the approach of continuous symmetry measures as developed by 
Zabrodsky, Pinsky and Avnir (2^j2^ might make it possible to bring together the surface and the volume symmetry 
properties. 
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